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We compute the graviton one loop contribution to a classical energy in a traversable wormhole 
background. Such a contribution is evaluated by means of a variational approach with Gaussian 
trial wave functionals. A zeta function regularization is involved to handle with divergences. A 
renormalization procedure is introduced and the finite one loop energy is considered as a self- 
consistent source for the traversable wormhole. 
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I. INTRODUCTION 

m . 
o ■ 

A wormhole can be represented by two asymptotically flat regions joined by a bridge. To exist, they must satisfy 
£NJ ' the Einstein field equations, like a black hole. While black holes are generally well accepted astrophysical objects, 
wormhole are yet to be discovered. One very simple and at the same time fundamental example of wormhole is 
represented by the Schwarzschild solution of the Einstein's field equations. One of the prerogatives of a wormhole 
is its ability to connect two distant points in space-time. In this amazing perspective, it is immediate to recognize 
i— ( ■ the possibility of traveling crossing wormholes as a short-cut in space and time. Unfortunately, although there is no 
direct experience, a Schwarzschild wormhole does not possess this property. It is for this reason that in a pioneering 
work Morris and ThornepJ and subsequently Morris, Thorne and YurtseverQ studied a class of wormholes termed 
11 traversable" . Unfortunately, the traversability is accompanied by unavoidable violations of null energy conditions, 
namely, the matter threading the wormhole's throat has to be "exotic" . Classical matter satisfies the usual energy 
conditions. Therefore, it is likely that wormholes must belong to the realm of semiclassical or perhaps a possible 
quantum theory of the gravitational field. Since a complete theory of quantum gravity has yet to come, it is important 
to approach this problem semiclassically. On this ground, Hochberg, Popov and Sushkov considered a self-consistent 
solution of the semiclassical Einstein equations corresponding to a Lorentzian wormhole coupled to a quantum scalar 
field||. On the other hand, Khusnutdinov and Sushkov fixed their attention to the computation of the ground state of 
a massive scalar field in a wormhole background. They tried to see if a self-consistent solution restricted to the energy 
component appears in this configuration Motivated by these works, we are interested to repeat the approach of 
Ref.Q but using gravitons instead of scalars. In particular, we wish to study the one loop contribution of the gravitons 
^ \ to the total energy. This is quite similar to compute the Casimir energy on a fixed background. It is known that, 
for different physical systems, Casimir energy is negative. This is exactly the feature that the exotic matter should 
possess. In particular, we conjecture that quantum fluctuations can support the traversability as effective source 
of the semiclassical Einstein's equations. The rest of the paper is structured as follows, in section [H] we define the 
effective Einstein equations, in section 11X11 we introduce the traversable wormhole metric, in section HV1 we give some 
of the basic rules to perform the functional integration and we define the Hamiltonian approximated up to second 
order, in section we study the spectrum of the spin- two operator acting on transverse traceless tensors, in section 
IVII we regularize and renormalize the one loop energy contribution and we speculate about self-consistency of the 
result. We summarize and conclude in section IVIII 
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II. THE EFFECTIVE EINSTEIN EQUATIONS 



We begin with a look at the classical Einstein equations 

Gfiu — K-T^v, (1) 

where T M „ is the stress-energy tensor, G^ v is the Einstein tensor and k = 8?tG. Consider a separation of the metric 
into a background part, g^, and a perturbation, h^ u , 

g^y = 9tlU + h^y. (2) 
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The Einstein tensor can also be divided into a part describing the curvature due to the background geometry 
and that due to the perturbation, 

G>„ {gap) = G>„ {(jap) + AG M „ (cjaf3, hap) , (3) 

where, in principle AG^ U (g a f3-,h a p) is a perturbation series in terms of h^ v . In the context of semiclassical gravity, 
Eq.Q becomes 

G>„ = k {T, v ) ren , (4) 

where (T liU ) ren is the renormalized expectation value of the stress-energy tensor operator of the quantized field. If the 
matter field source is absent, nothing prevents us from defining an effective stress-energy tensor for the fluctuations 
as 1 

(T^y en = -- (AGpu (g a0 , h aP )) ren . (5) 

From this point of view, the equation governing quantum fluctuations behaves as a backreaction equation. If we fix 
our attention to the energy component of the Einstein field equations, we need to introduce a time-like unit vector 
such that u ■ u = — 1. Then the semi-classical Einstein's equations Q projected on the constant time hypersurface E 
become 

G>„ {g aP ) u»u v = k (T^u») ren = - (AG M „ (g a/3 , h a p) u»u v ) ren . (6) 
To further proceed, it is convenient to consider the associated tensor density and integrate over E. This leads to 2 



where 



^ / d 3 xV^G^ (g af3 ) uV = - [ d 3 xH i0) = ~ ( d 3 xy/^ (AG^ (g af3 , h af) ) u^u") ren , (7) 



v 9 K 

is the background field super-hamiltonian and G^i is the DeWitt super metric. Thus the fluctuations in the Einstein 
tensor are, in this context, the fluctuations of the hamiltonian. To compute the expectation value of the perturbed 
Einstein tensor in the transverse-traceless sector, we use a variational procedure with gaussian wave functionals. In 
practice, the right hand side of Eq.JTJ) will be obtained by expanding 

(tf|ff E |tf\ h H<®+H<P+H<P 



(9) 



and retaining only quantum fluctuations contributing to the effective stress energy tensor. H^r represents the hamil- 
tonian approximated to the i th order in and $ is a trial wave functional of the gaussian form. Then Eq.JjJ 
becomes 

4 13 +4 2) + 



= / d 3 xH^ = - ' * 



(10) 



The chosen background to compute the quantity contained in Eq.Q will be that of a traversable wormholc. 



1 Note that our approach is very close to the gravitational geon considered by Anderson and Brill|(j. The relevant difference is in the 
averaging procedure. 

2 Details on sign conventions and decomposition of the Einstein tensor can be found in ApeendixlAl 
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III. THE TRAVERSABLE WORMHOLE METRIC 

In Schwarzschild coordinates, the traversable wormhole metric can be cast into the form 

dr 2 

ds 2 = - exp (-20 (r)) dt 2 + ^ + r 2 [d9 2 + sin 2 6dip 2 ] . (11) 

where <f> (r) is called the redshift function, while b (r) is called the shape function. Proper radial distance is related to 
the shape function by 



f r dr' 

l{r) = ± -J^vT. 



(12) 

where the plus (minus) sign is related to the upper (lower) part of the wormhole or universe. Two coordinate patches 
are required, each one covering the range [r t , +oo). Each patch covers one universe, and the two patches join at rt, 
the throat of the wormhole defined by 

r t = min{r(0}. (13) 
Instead of Ea. ((TT|l . we consider the special case where cb = and b (r) = r 2 /r. Then, we get 

ds 2 = -dt 2 + + r 2 [dO 2 + sin 2 0d V 2 } (14) 

1 - £ 

and with the help of Ea. (|12|l . we obtain 

dr 1 i 

r 2 = I 2 + r\ (15) 



f r dr' r 



and Ea. 1)14(1 becomes 

ds 2 = -dt 2 + dl 2 + (r 2 + I 2 ) [dO 2 + sin 2 9dtp 2 ] . (16) 

The new coordinate I covers the range — oo < I < +oo. The constant time hypersurface E is an Einstein-Rosen bridge 
with wormhole topology S 2 x R 1 . The Einstein- Rosen bridge defines a bifurcation surface dividing £ in two parts 
denoted by S + and To concretely compute Ea. l|lU|) . we consider on the slice E deviations from the wormhole 
metric of the type 

9ij =9ij + fHj, (17) 
where §ij is extracted from the line element l(16|) whose form becomes 

ds 2 = -dt 2 + t i, ,(l.i-'d.i- ! . (18) 

IV. ENERGY DENSITY CALCULATION IN SCHRODINGER REPRESENTATION 

In order to compute the quantity 

' d 3 x V^ (AGV (g a0 , h a0 ) u^) ren , (19) 



we consider the right hand side of Ea. (|10|l . Since is linear in and h, the corresponding gaussian integral 

disappears and since 

y G llv {g aP )u^u v = -H, (20) 
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it is clear that the hamiltonian expansion in Eq. i|10|) does not coincide with the averaged expanded Einstein tensor of 
Eci. (|19fl because Eq.J^UJ) involves a tensor density. Therefore, the correct setting is 



d 3 x^(&G^(g a p,h al3 )u^) ren = / d 3 x^g 



(*l*> 



"(2) . 



where \/ 3 g is the second order expanded tensor density weight. Following the same procedure of Refs. 
potential part of the right hand side of Ea. l|21|) becomes 



/ d 3 x^ 



The term 



-hAh + -h ll Ah H - -h ij ViVM + -h\7i\7ih u - -h ij R ia h<> + -hR t ,h lj + -Rh ij h u - -hRh 

A A 9 J o 9 J o J A L J A 



/ d 3 Xy/^ 



-Rh i] h^ - -hRh 
4 3 4 



(21) 
, the 

(22) 
(23) 



makes the difference between the hamiltonian expansion and the Einstein tensor expansion. To explicitly make 
calculations, we need an orthogonal decomposition for both and hij to disentangle gauge modes from physical 
deformations. We define the inner product 



(h,k) := / ^gG ijkl h ij (x)k k i(x)d 3 x, 



(24) 



by means of the inverse WDW metric Gijki , to have a metric on the space of deformations, i.e. a quadratic form on 
the tangent space at h, with 



Qijkl _ (g*k g jl + gilgjk _ 2^"). 

The inverse metric is defined on co-tangent space and it assumes the form 



(25) 



(P. l) ■= / \ gGijkip ij (x) Q kl (x) d 3 x, 



so that 



(26) 



(27) 



Note that in this scheme the "inverse metric" is actually the WDW metric defined on phase space. The desired 
decomposition on the tangent space of 3-metric deformations |0 ^3 i s: 



hij = + i L Oij + h ij 

where the operator L maps £j into symmetric tracefree tensors 



(28) 



= + V 3 -& - -9ij (V • . 



(29) 



Thus the inner product between three-geometries becomes 



(h, h) := / y /gG ijkl h ij (x)h kl (x)d 3 x 



V9 



-Ih^ + myHLOij + h^htj 



(30) 
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With the orthogonal decomposition in hand we can define the trial wave functional 



V{hij (a?)} = Afexp 



1 

4Zl 



-1 l \ Trace 



(31) 



where J\f is a normalization factor. We are interested in perturbations of the physical degrees of freedom. Thus we 
fix our attention only to the TT tensor sector reducing therefore the previous form into 



*{h ij (£)}=Mexp{-^(hK- 1 h)^ i 



(32) 



Therefore to calculate the energy density, we need to know the action of some basic operators on VP [hij]. The action 
of the operator hij on \9) = 9 [hij] is realized by 



hij (x) |*) = hij (It) ^ {h^}. 



(33) 



The action of the operator mj on j^), in general, is 



TTij (x) \9) 



-9 {hij}. 



Shij (~x) ~ L ' UJ 
The inner product is defined by the functional integration: 

(* i | *a) = / Phj] 9{ {hij}9 2 {h kl } , 



(34) 



(35) 



and by applying previous functional integration rules, we obtain the expression of the one-loop-likc Hamiltonian form 
for TT (traceless and transverse) deformations 



d 3 x^G 



ijkl 



[1G-kG)K- ix {x,x) ijkl + 



The propagator (x,x) iak[ comes from a functional integration and it can be represented as 



ikl 



E 



ftj r)± (^)4l u (^) 

2A(r) 



(36) 



(37) 



where h^ ± (~x) are the eigenfunctions of A2, whose eigenvalues will be denoted with to 1 (r). r denotes a complete 
set of indices and A (r) are a set of variational parameters to be determined by the minimization of Eq. Ij36[l . The 
expectation value of H 1 - is easily obtained by inserting the form of the propagator into Eq. H36JI 



(16ttG) A, (r) + 



"i 2 (r) 



167rG) A, (t) 



(38) 



By minimizing with respect to the variational function A^ (r) , we obtain the total one loop energy for TT tensors 



(39) 



The above expression makes sense only for u>f (r) > 0, i = 1, 2. The meaning of ujf will be clarified in the next section. 
Coming back to Ea. (|10|l . we observe that the value of the wormhole energy on the chosen background is 



and to one loop the self-consistent equation for TT tensors becomes 

2G = ~ E ■ 



(40) 



(41) 



Note that the self-consistency on the hamiltonian as a reversed sign with respect to the energy component of the 
Einstein field equations. This means that an eventual stable point for the hamiltonian is an unstable point for the 
effective energy momentum tensor and vice versa. 
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V. THE TRANSVERSE TRACELESS (TT) SPIN 2 OPERATOR FOR THE TRAVERSABLE 
WORMHOLE AND THE W.K.B. APPROXIMATION 

In this section, we evaluate the one loop energy expressed by Ea. (|39|l . To this purpose, we begin with the operator 
describing gravitons propagating on the background . The Lichnerowicz operator in this particular metric is 
defined by 



A 2 h^ : = - (A T h^)l + 2 (Rh^ - R (fc") , 



(42) 



where the transverse-traceless (TT) tensor for the quantum fluctuation is obtained by the following decomposition 

(43) 



I-*. 1.,-. ,,3V 1 



hl=hl--6lh+-6lh=(h T y i + -Slh. 
This implies that (h T Y. = 0. The transversality condition is applied on (h T Y and becomes Vj (h T Y- = 0. Thus 



faint = - a s {h TT )\ + £ 1 1 



b(r) 



where A 5 is the scalar curved Laplacian, whose form is 



b(r) 



4r - b' (r) r - 36 (r) 
2V 2 



dr r 2 



and R°j is the mixed Ricci tensor whose components are: 



pa 



b'(r) b(r) b'(r) b{r) b' (r) 6 (r) 



The scalar curvature is 



2r 2 



R = R{8) 



2r 3 ' 2r 2 



,6'(r) 



2r 3 



We are therefore led to study the following eigenvalue equation 



(44) 



(45) 



(46) 



(47) 



(48) 



where uj 2 is the eigenvalue of the corresponding equation. In doing so, we follow Rcggc and Wheeler in analyzing 
the equation as modes of definite frequency, angular momentum and parity [l^. In particular, our choice for the 
three-dimensional gravitational perturbation is represented by its even-parity form 



diag 



H(r)[l 



b(r) 



, r 2 K (r) , r 2 sin 2 <&L (r) 



Yi m (■&, 0) . 



(49) 



with 



H(r) = h\ (r)-U(r) 
K{r) = hl (r)-U(r) 
L{r)=hl (r)-lh(r) 



(50) 



From the transversality condition we obtain h 2 , (r) = (r) . Then K (r) = L (r) . For a generic value of the angular 
momentum L, representation (|49|l joined to Ea. l|44(l lead to the following system of PDE's 



-A, + 2 



(-A,+2( 



fc'(r-) _ b(r) _ b'(r) 



b'(r) , b(r) _ b'(r] 



2r 2 



2r 3 



)H{r)=u 2 l H{r) 
)K(r)=<4jK(r) 



(51) 



where A; is 



A, = 1 



Defining reduced fields 



b(r)\ d 2 f4r-b'(r)r-3b(r)\ d 1(1 + 1) 6 / b(r) 



r J dr 2 \ 



2r 2 



dr 



r r 

and passing to the proper geodesic distance from the throat of the bridge defined by 

dr 



the system ()51f) becomes 



dx = ± 



f 



b(r) 



(52) 



(53) 



(54) 



with 



■i? + V l {x) h{x)^Lo 2 J 1 {x) 



■£ + V 2 (x)]f 2 (x)=u>y 2 (x) 



V 1 (x) = ^ + U 1 (x) 
V 2 (x) = ^ + U 2 (x) 



and 




(55) 



(56) 



(57) 



In order to use the WKB approximation, we define two r-dependent radial wave numbers k\ (x, I, u>i t ni) and 
k 2 (x,l,oj 2>n i) 



k\ (ar, I, wi.ni) = u\ nl - $X% ~ Ul ( x ) 

The number of modes with frequency less than Wi, i = 1, 2, is given approximately by 

Ufa) = J Vi{l,Ui)(2l + l)dl, 
where (I, LOi), i = 1,2 is the number of nodes in the mode with (I, oji), such that 

Vi (l,U!i) = — [ dxJ fcf (X,l,t0i). 

n J-oc v 



(58) 



(59) 



(60) 



Here it is understood that the integration with respect to x and I is taken over those values which satisfy kf (x,l,0Ji) > 
0, i = 1, 2. Thus the total one loop energy for TT tensors is given by 



E 



TT 
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where 

' Pi 



+ 00 

dx (x 2 + r 2 t ) [p x + p 2 ] , (61) 



(62) 



VI. ONE LOOP ENERGY REGULARIZATION AND RENORMALIZATION 



In this section, we proceed to evaluate the one loop energy. The method is equivalent to the scattering phase shift 
method and to the same method used to compute the entropy in the brick wall model. We use the zcta function 
regularization method to compute the energy densities p\ and p 2 - Note that this procedure is completely equivalent to 
the subtraction procedure of the Casimir energy computation where zero point energy (ZPE) in different backgrounds 
with the same asymptotic properties is involved. To this purpose, we introduce the additional mass parameter p, in 
order to restore the correct dimension for the regularized quantities. Such an arbitrary mass scale emerges unavoidably 
in any regularization schemes. Then we have 



1 f + °° 

Pi ( £ ) = / , du; 

47r JJutM 



V^) (uj 2 - Ui (x)) e 



(63) 



If one of the functions Ui (x) is negative, then the integration has to be meant in the range where uj 2 + Ui (x) > 0. In 
both cases the result of the integration is 3 



64tt 2 



e 



P 



Ui(x) 



2 In 2 - - 
2 



(64) 



where the absolute value has been inserted to take account of the possible change of sign. Then the total regularized 
one loop energy is 



E TT (r u e;p) 



dx (x 2 + r 2 ) [( Pl (e) + p 2 (e))] . 



The result of the integration over the x coordinate leads to the following expression 



E TT (r t ,e;u) = — 
1 ' ,w 16tt 



a b c 

In (r t p) 

er t r t r t 



(65) 



(66) 



where a = 23.35740218, b = 276.6026775 and c = 212.0575042. Then the self consistent equation @TJ can be written 
in the form 



2G 



1 

16tt 



a 

en 



b 



In (r t p) 



(67) 



In order to deal with finite quantities, we renormalize the divergent energy by absorbing the singularity in the classical 
quantity. In particular, we rc-dcfinc the bare classical constant G 



1 

G 



1 

Gq 



e87r 2 r t 2 ' 

Therefore, the remaining finite value for the effective equation (|67() reads 



2G 



1 

167T 



b c 1 , \ 
1 In (r t p) 

n n 



(68) 



(69) 



3 Details of the calculation can be found in the Appendix IH1 
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This quantity depends on the arbitrary mass scale m- It is appropriate to use the renormalization group equation to 
eliminate such a dependence. To this aim, we impose that|l4| 



d 



irr t 



1 



namely 



d/i \ 2Go (m) J dfJ, ( 167T 

nr t OGq 1 (/i) c 



& c i i \ 
1 In [rtli) 

n n 



2~ M ~ d^i 



0. 



(70) 



(71) 



Solving it we find that the renormalized constant Go should be treated as a running one in the sense that it varies 
provided that the scale \i is changing 



1 



or G (m) 



G (mo 



K = 



G (fi) G (mo) ' " V/V "' " uv ^ l + G (Mo)^ln(^)' " 8 ^ r * 2 ' 
where Mo is the normalization point. We substitute Ea. (|72(l into Ea. H69() to find 



(72) 



2G (mo) 16tt 



4- + 4" m ( r *Mo) 



(73) 



where we have divided by r t . In order to have only one solution 4 , we find the extremum of the r.h.s. of Eq.J73J) and 
we get 



2b 



2c 



In (r t Ho) 



1 (c~2b 
n = — exp — — 
Mo V 2c 



and 



Go (/xo) 167T 
With the help of Eqs.dJ and f75|l. Eq.O becomes 

1 - 1 ■ C ln^ 



CMn /" c-2b 
— exp 



G (m) Go (mo) ' 87r 2 rf 



Mo 



cun /" c — 2b 
T7T^ exp 



1 + 2 In 



Mo 



(74) 



(75) 



1 



l + 21n[ 

Mo 



Go (Mo) 

It is straightforward to see that we have a constraint on m/Mo- Indeed we have to choose 

r 



M > Mo exp 



.6065306597mo, 



otherwise Go (m) becomes negative. We have now two possibilities: 

1. we identify Go (mo) with the squared Planck length. As a consequence we obtain 



f t = —y/cGo (mo) = 1.158822606Z P 

47T 



Mo = .3860531213m p . 



(76) 



(77) 



(78) 



4 Note that in the paper of Khusnutdinov and Sushkov 4] , to find only one solution, the minimum of the ground state of the quantized 
scalar field has been set equal to the classical energy. In our case, we have no external fields on a given background. This means that it 
is not possible to find a minimum of the one loop gravitons, in analogy with Ref.0- Moreover the renormalization procedure in Ref.0 
is completely independent by the classical term, while in our case it is not. Indeed, thanks to the self-consistent equation 1411 . we can 
rcnormalize the divergent term. 
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2. We identify /Uo with the Planck scale and we get 

f t = .4473670842^ ^/G (/i ) = .3860531213Z P . (79) 

If we compare our result with that obtained by Khusnutdinov and SushkovQ, we see that in their work they 
find that the wormhole radius is 

r w ~ .0U11 P (80) 

and its mass 

m w ~ 11.35m p (81) 

is trans-planckian. We recall that these values have been obtained by looking at a massive scalar field quantized 
in a wormhole background playing the role of the exotic matter. In our case, we suppose that the graviton 
quantum fluctuations play the role of the exotic matter and even if we fix the renormalization point at the 
Planck scale as in Ref.Jj, we find a radius f t > r w . 



VII. SUMMARY AND CONCLUSIONS 



Motivated by the works of Anderson and Brill, in searching for gravitational geons|6j and of Khusnutdinov and 
SushkovQ, we have considered the possible existence of a self-consistent solution of the semiclassical Einstein's 
equations in a traversable wormhole background. In particular, we have fixed our attention to the graviton quantum 
fluctuations around such a background. The fluctuations, contained in the perturbed Einstein tensor, play the role 
of the "exotic" matter considered in Ref.[4|. A variational approach with the help of gaussian trial wave functionals 
has been used to compute the one loop term. The adopted procedure is quite close to that of Anderson and Brill 
except for the averaging process which, in our case, involves a variational calculation. To handle the divergencies we 
have used a zeta function regularization which, in this context, is formally equivalent to a Casimir energy subtraction 
procedure. A renormalization of the Newton's constant G has been performed to absorb the ultraviolet divergencies. 
To avoid dependences on the renormalization scale a renormalization group equation has been computed. While the 
renormalization process is not new in the context of semiclassical Einstein field equations, to our knowledge it is the 
use of the renormalization group equation that seems to be unknown, especially concerning self-consistent solutions 
and traversability. The interesting point is the obtained consistency with a wormhole radius f t greater than radius 
r w of Ea. (|80|l . It is likely that this result be correlated with the graviton themselves. Despite of this, the obtained 
11 traversability" has to be regarded as in "principle" rather than in "practice" because of the obtained wormhole 
radius size. Nevertheless, the calculation is far to be complete because only gravitons have been taken into account 
and the traversable wormhole metric considered is restricted to Ea. (|14fl . A possible improvement of our evaluation 
could come from including a charge[l5j or by looking at a Schwarzschild-like metric, but even with a different choice 
of the wave functionals entering in the variational approach. 



APPENDIX A: EINSTEIN EQUATIONS AND THE HAMILTONIAN 

Let us consider the Einstein equations 

G^u = R^ u - \^9p.vR = k8ttGT^. (Al) 

Rix V is the Ricci tensor and R is the scalar curvature. If it M is a time-like unit vector such that g^u^u 1 ' = — l,then 
the Einstein tensor G^y becomes 

G^vTu" = R^vTu" - X -g^u v R = R^u" + ^R. (A2) 

By means of the Gauss-Codazzi equations^, 

R = i? (3) ± 2R flu u>*u v tK 2 ± K^K^, (A3) 

where is the extrinsic curvature and R^ is the three dimensional scalar curvature. For a time-like vector, we 
take the lower sign and Ea. (jA2|) becomes 

G^uW = \ (i? (3) + K 2 - K^K^ . (A4) 
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If the conjugate momentum is defined by 



/(3) 



2k 
then 

2 



(Kg** - K» v ) , (A5) 



K 2 -K^K^=(^=) l^-ir^) (A6) 



/(3) <? 



and 



^„ ., „ \/^Q ^ 2k fn 



G^u v = ^i?( 3 ) + ^ - 7r"%„ - , (A7) 



2k ^ ~ 2k " ' ^(3)^ V 2 
namely Eq.JSJ with the reversed sign. 

APPENDIX B: THE ZETA FUNCTION REGULARIZATION 

In this appendix, we report details on computation leading to expression (|63() . We begin with the following integral 

I + =^f+ 00 dcj £ r 

P{e)={ , (Bl) 

I_ = p** f+°° fa ^ r 

•"■> (,.,2_TT(~\\e-n 



with U (x) > 0. 



1. I + computation 



p ( e ) = M 2e [/ ( x ) 2 ~ £ / dt r = t:^U {xf~ e B [ s - 2 



If we define t = u)/y/U (x), the integral 7 + in Ea. ljBip becomes 

t 2 _ 

where we have used the following identities involving the beta function 

ft _|_ QQ -^2X~~\ 

B(x,y)=2 dt Rex>0,Rey>0 (B3) 

Jo ft 2 + 1 y 



related to the gamma function by means of 



r (x + y) 

Taking into account the following relations for the T-function 



B^y) = T -^4. (B4) 



r(e - 2) = e(e -l) (e l2) ' H^^^f' ^ 



and the expansion for small e 



r ( i + ) = i - + ( > ( 2 i r ( £ + = r (l\ - eT (I ) ( 7 + 2 in 2) + o ( ■-) 



12 



x £ = 1 + ehxx + O (e 2 ) , 



where 7 is the Euler's constant, we find 



U 2 (x) 
16 



- + In 



[i 1 \ 1 



(B6) 
(B7) 



2. /_ computation 



If we define t = oj/ y/U (x), the integral I- in Ea. ljBlfl becomes 

t 2 1 



P (e)=^U{x) 



\^U(x) 2 -^ 



-\-oo 



. T = L^U{x) 2 - e B[e-2r--e 
(£ 2 -l) e ~ 5 2 V 2 



where we have used the following identity involving the beta function 



p V p 



+oc 



dtt^ 1 (t p - 1)" 



p > 0, Rc v > 0, Re ^ < p — p Re v 



and the reflection formula 

r (z) r (1 - z) = -zv (-z) r (z) 

From the first of Eas. (|B5|> and from the expansion for small e 



r (|- £ ) = r (|) (l-e(-7-21n2 + 2)) + 0( £ 2 ) 



1 + elnx + O (e 2 ) , 



we find 



P(e) = - 



uHx)_ 

16 



- + In 



U x) 



2 In 2 - - 
2 



(B8) 

(B9) 
(BIO) 



(BH) 
(B12) 
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